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Abstract 

Dense granular media have a compaction between the random loose 
and random close packings. For these dense media the concept of a gran- 
ular pressure depending on compaction is not unanimously accepted be- 
cause they are often in a "frozen" state which prevents them to explore 
all their possible microstates, a necessary condition for defining a pressure 
and a compressibility unambiguously. While periodic tapping or cyclic 
fluidization have already being used for that exploration, we here suggest 
that a succession of flowing states with velocities slowly decreasing down 
to zero can also be used for that purpose. And we propose to deduce the 
pressure in dense and flowing granular media from experiments measuring 
the thickness of the granular layer that remains on a rough incline just 
after the flow has stopped. 



1 Introduction 

The existence of a pressure in granular media at rest is the simplest way to 
represent their stiffness. When the concentration of grains is above the random 
close packing, the granular medium acts as a poro-elastic solid and its pressure 
is a function not only of the compaction but also of the elastic constants of 
the material the grains are made of. When the granular medium has a smaller 
compaction, in the range between the random loose and random close packings, 
the grains can be considered as rigid and the expression of the granular pres- 
sure is far less evident. The difficulty comes from the glassy behaviour which 
makes it quite usual to find dense granular media in a frozen state concern- 
ing their compaction. Hence the feeling that the granular pressure is largely 
dependent on the way the medium was prepared. However, experiments have 
been conducted which aim at allowing the granular medium to reach a steady 
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and quasi-equilibrium state concerning its compaction. These experiments re- 
lied on regular tappings or cyclic fluidization favouring the exploration of many 
microstates (see e.g. £Q and [2]). We propose that a systematic exploration of 
the microstates can also be achieved starting from a flowing granular medium, 
and slowly reducing its velocity down to rest. As a consequence, we suggest 
that some of the experiments on rough plates (inclined with angle 9) which led 
to define the thickness h s top{9) which remains after the flow has stopped, can 
also be used to infer the relation between granular pressure and compaction. 
While the pressure in poro-elastic media originates from elastic forces and gran- 
ular deformations, the pressure in dense and flowing granular media (those with 
a compaction between the random loose and random close packings) results 
mainly from the many possible random spatial configurations of the grains. We 
will present the main features of this "randomness" pressure as well as those 
of the more classical elastic pressure. Then some general expressions for the 
granular pression will be proposed. Finally two particular expressions will be 
deduced from experimental data concerning h stop - 

2 The Configuration Pressure and the Disorder 
Pressure 

Consider a large volume containing many rigid spheres with a high enough vol- 
ume fraction <j> for a contact network to invade the whole volume. Let Q(<p)d<p 
be the number of different spatial configurations of these spheres in the range 
between <f) and (f> + dcf>. To belong to f2($), a configuration must display a large 
enough number of contacts, but with zero forces at the contact points. The 
density of micro-states f2(0) is thus a purely geometric concept. If instead of 
rigid spheres we were considering soft ones, we would say we are counting the 
number of configurations with zero energy ("incipient" contacts), yet able to 
resist an infinitesimal external pressure load. This density of states presumably 
vanishes below a minimum compaction cj) m in — 0.40 (the gel threshold) and 
above the maximum compaction 4> max — 0.74 (the most compact crystalline 
configuration). At some intermediate volume fraction cj> m the density of mi- 
crostates displays a large maximum value. This intermediate compaction with 
the maximum number of microstates happens to be the loosest random packing. 
According to Onoda and Liniger :'> its value for spherical grains is (f> m = 0.555. 
Introduce now a number P which represents a non-dimensional measure of the 
configuration pressure and define the partition function 

Z{P) = e p ^ il(<f>) d<f>. 

The mean volume fraction is then related to the configuration pressure in the 
form: <f> = -jp Log Z. This relation can be used to obtain both P(cf>) and the 
variance of the density fluctuations < ((f) — <fi) 2 >= d(f)/dP. Another useful 
quantity is the configuration entropy S — Log Z — P<p . This configuration 
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entropy is a function of the mean volume fraction and the non-dimensional 
configuration pressure is nothing but P = —dS/dcj>. Many works starting 
from Kanatani [3] and revived by Edwards and Oakeshott [Sj strived to find an 
explicit form for the partition function or for the configuration entropy. The 
general trends of their results are the followings : for a vanishing pressure, the 
mean volume fraction is <j) m while for an infinite pressure, the volume fraction is 
4>max and the compressibility d(j)/dP vanishes. A very simple expression meeting 
these conditions is 

configuration 



p con figuration ^ p ^ Log 

4*max 4* 

Close to the maximum volume fraction the configuration pressure has an expres- 
sion quite similar to the pressure deduced from the entropy of the lattice-gas 
model (sec e.g. 6 ). Note that the configuration pressure stems from the total 
number of different configurations, including both random and crystalline ones. 
It is also possible to select random configurations only. In this case one intro- 
duces the random close compaction <f>M above which all configurations display 
some cristalline order. For spheres, it is generally admitted that </>m = 0.635 
and a possible expression for the disorder pressure is 

^disorder ^ £ g < ^ >M ~ < ^ >m (1) 

<t>M — 4> 

The gradient of the disorder pressure acts as a diffusion force which pushes the 
grains towards lower compactions, those with a larger number of microstates. 
For the grains to have a chance to explore all microstates with equal probability, 
the best solution is a steady flow. This is why the concept of a disorder granular 
pressure is more pertinent for dynamic situations and that it must be handled 
with care in static ones. The disorder pressure confers the granular medium a 
compressibility (variance of the compaction fluctuations) which decreases with 
the compaction. When taking the above expression for granted the compress- 
ibility is proportional to <pM — 0, in agreement with the experimental results 
of Nowak et al. but not with those of Schroter et al. 2 which display a 
minimum of the compressibility for a compaction between <j> m and 4>m- The 
main features of p dtsorder are drawn schematically in Fig.l. 



3 The Elastic Pressure 

The configuration or disorder pressure is a purely geometric (or entropic) con- 
cept. What matters is a minimum coordination number between particles but 
the forces at the contact points are of no concern. It is then obvious that the 
configuration pressure is not the whole story. There is a second (and more in- 
tuitive) source of granular pressure which witnesses to the elastic stiffness of 
the granular medium and for which the magnitude of the contact forces is of 
utmost importance. That second contribution is the elastic pressure p elastlc . 
According to Aharonov and Sparks [B] and to O'Hern, Silbert, Liu and Nagel 
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P| (see also Head and Doi ^Uj for a first principles approach) p elastlc vanishes 
below a volume fraction (ft c which depends on the friction between grains and 
is sligthly smaller than the random close packing (ft>M- Above (ft c it was found 
that p elastlc (0 — c ) a with a = 1 for Hookean contact forces and a — 3/2 for 
Hertzian ones. What is not apparent in the above expression is that the bulk 
modulus dp elasUc / d(ft is usually orders of magnitude larger than dp dlsorder / d(ft. 
When plotted on the same graph as p dlsorder ^ it is as if p elasUc was zero for 
(ft < <p c and infinite for (ft > <p c> as sketched in Fig.l where the difference be- 
tween <ft c and 4>m was exagerated for clarity. 



4 The Rate-Independent Granular Pressure 

The pressure in a moving granular medium is generally made of two parts : 
a rate-dependent part which represents dilatancy effects in shear flows and a 
rate-independent part (hereafter noted p) representing the physical phenomena 
involved in both p elastlc and p dlsorder _ We are interested in the latter rate- 
independent contribution p. It seems clear that the very high stiffness for (ft ~ (ft c 
is a feature of p elasttc that must be shared by the total granular pressure. Much 
less obvious is the behaviour of p for volume fractions slightly above (ft m . Should 
we trust the linear behaviour displayed by the above expression for p dlsorder ? 
In fact, since our final aim is a comparison with experimental results, it seems 
wise to introduce a general expression like p ~ {(ft — (ft m ) n /(4>m — <ft) m where 
the positive exponents m and n are left undetermined and we neglected the 
small difference between <ft c and 4>m- Note that this expression for the rate- 
independent granular pressure holds in the very small range (ft m < (ft < <j>M 
specific of dense granular media, and that the granular pressure vanishes for 
(ft = 4> m and diverges for (ft — <ft>M. In what follows we scale that granular- 
pressure with p p gD where p p is the mass per unit volume of the grains and 
g is the acceleration of gravity. That scaling seems obvious when gravity is 
responsible for the confinment of the grains but we checked that the same 
scaling is also convenient when the confinment is due to an external pressure 
applied on the granular material. The physical meaning of that special scaling 
is clear : the loads exerted on a granular medium with compaction in the range 
4> rn < (ft < 4>m are much smaller than those leading to poro-elastic media with 
(ft > (f>M- The latter ones involve elastic deformations of the grains while grains 
with a smaller compaction can be considered as rigid. Consequently we can not 
accept any scaling involving elastic constants and we are left with p p gD as the 
only relevant order of magnitude for stresses. Defining the relative compaction 
(p as 

4>M - <ftm 

we will test the potentialities of two expressions for the rate-independent gran- 
ular pressure in the range < ip < 1 : a rather general one with exponents 
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m > and n > 

P = P Pp9 D (1 _ , (3) 
and a special one identical to p dlsorder a s given in JTJ 

P = Po P P gD Log — - (4) 

That latter expression was already adopted in previous works [121 ll.'ij and it 
suggests that the elastic pressure plays a negligible role except for <fi ~ </>jv/. We 
now consider an experiment likely to discriminate between lyj and and able 
to give an order of magnitude for Pp. 



5 The Maximum Thickness of a Granular Layer 
that Stops on a Rough Incline 

Consider a layer of granular material flowing down a rough inclined plate. Upon 
gently reducing the inclination, the layer ultimately stops at some angle 9, with 
a thickness h. Since the flow velocity was slowly reduced to zero, the granular 
layer had time to explore a lot (if not all) of the microstates involved in p disorder _ 
And since the flow was slow, the rate-dependent part of the granular pressure 
was already negligible before the layer jams. It is thus likely that the peculiar 
jammed state which the layer arrives at is described by the rate-independent 
pressure p defined above. The mechanical equilibrium of the freshly jammed 
layer is thus expressed by 

dp 

= — — — h 4>p p gcos6, and tan(9 = min[fi(z)], 

where the z— axis is orthogonal to the free-surface of the layer and points down- 
wards while p is the friction coefficient which is possibly non-uniform over the 
layer thickness. Substituting expression (|3J), one deduces that the reduced com- 
paction profile ip(z) increases from zero close to the free-surface up to values 
close to one at a distance of order L with -fe = , p " „ . Close to the free-surface 

JJ (pMCOSti 

that is to say for < z <C L, and far from the free-surface [L <C z) the profile 
behaves asymptotically like <p(z) — (^f;)" an d v{ z ) = 1 ~ ("?)"'■ Concerning 
the particular pressure the whole profile is exponential-like and given by 

A similar " Fermi-Dirac" profile was already observed in experiments |14j and 
was also obtained in simulations of a frustrated lattice-gas model |15| . 
Knowing the compaction profile, let us now focus on the layer thickness h. Da 
Cruz has deduced from numerical simulations a very important (and a bit 
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counter-intuitive) result : the macroscopic friction coefficient decreases almost 
linearly with the compaction and can be written as 



fi = tan9 max - {tanOmax - tan6 min ) (p . (6) 

Since the compaction increases with the distance from the free-surface, the 
smallest value of fi happens very close to the rough plate, that is to say for z ~ h. 
Just after jamming, we thus have tan9 max — (tan9 

max tO/TlOfjilYi 

tan#. Taking the compaction profiles (??) into account, one deduces that when 
9 slightly exceeds 6 min the layer is very thick and 



h s top Pq I tan 9 rnax — tan I 



D 6 Mcos0 m i n \ tan 9 tan^^^ 
while when 9 comes close to 9 max the layer is thin and 
h s top _ Pq ( tan 9 max - tan ( 



D <t> m cos6 max V tan 6 max - tan 6, 



(7) 



(8) 



For the special case (@J, one obtains a rather simple expression that holds in 
the whole range 9 min <9< 9 max 



hstop Po T 
= Log 

D 4>mcos9 



1 4>m tan 9 max - tan t 
6 m tan 9 - tan 9 rnir 



(9) 



6 Experimental results 

Systematic measurements of the layer thickness were initiated by Pouliquen |17j 
who distinguished between the thickness h s t a rt(9) for an initially static layer and 
h s top{9) for an initially flowing layer. The experimental results were fitted by 
two different expressions 

hstart,sto P (9) „ T tan 6*2 — tan#i h s tart,sto P {9) tan 6*2 — tan# 

[ — = B Loq or — = B . 

D L tan# - tan 0i D tan6>-tan#i 

(10) 

where D is the grain size while B, 9\ and 9^ are constants, different for h s top 
and hstart- Assuming that 9i — 9 max and 9\ = 9 min , these data fittings sug- 
gest two possible expressions for the granular pressure : either expression Q 
or expression @ with n — 1 and m = 1. However, most of the experimental 
values obtained for B, 9i and #2 must be considered as non representative of 
the bulk granular pressure because they were found to depend on the roughness 
of the inclined plate. In fact, the main difficulty with the interpretation of ex- 
perimental results is to estimate the effective friction coefficient associated with 
this roughness. Our theoretical predictions were based on one main assump- 
tion : the minimum of effective friction occurs somewhere close to the incline, 
at the limit between the bulk and a thin boundary layer strongly influenced 
by the plate roughness. This supposes the effective friction in the boundary 
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layer to be much larger than its value in the bulk. In other words, to deduce 
the bulk granular pressure from experiments on inclines, we are led to exclude 
those experiments performed with relatively smooth plates, and more generally 
those for which the curves h sta rt,stop(@) are strongly modified upon changing 
the plate roughness. And concerning those with a high enough roughness, we 
must exclude some boundary layer of thickness S and consider h — S as the rel- 
evant thickness for bulk behaviour. Accordingly, we were led to discard all the 
experiments performed with glass beads because the friction generated by the 
beads glued on the incline is only slightly larger than the friction in the bulk. 
But we considered as significative the experiments with sand flowing on carpets 
of various roughnesses. And for these experiments with sand we discarded a 
boundary layer with thickness estimated to 6 ~ AD. Because of the scarcity 
of data, we could not fully discriminate between expression (0J and expression 
with m = 1 and n = 1. We could fit the experimental results for sand on 
carpets (see Fig. 2) with 

[a = 0.73 — Q.l6(p and p~ 5.3 p„gD Log , (11) 

as well as with 

H = 0.73 - 0.20</? and p ~ 5.2 Pp gD . (12) 

7 Conclusions 

The concept of a rate-independent granular pressure was proposed long ago 
for poro-elastic media We suggest this concept can be extended to dense 
and flowing granular media for which this compaction-dependent pressure is 
the consequence of two distinct physical phenomena, of entropic and mechanis- 
tic nature respectively. We proposed to deduce the bulk granular pressure from 
experiments on rough inclines and particularly those giving h stop (8) . It happens 
that most of the experimental results concerning this thickness depend more on 
the roughness of the incline than on the bulk behaviour of the granular medium. 
However, the scarce results representative of the bulk behaviour can be inter- 
preted with a granular pressure and an effective friction depending both on the 
compaction. Concerning sand, a rather satisfactory fit with the experimental 
results could be obtained with either expression (|llf> or expression (|12f) . wit- 
nessing to the predominance of entropic effects in almost the whole compaction 
range <f> m < 4> < 4>m , while elastic effects come into play when approaching the 
random close packing cf>M only- 
One may wonder why numerical simulations of flows down rough inclines (see 
e.g. predict flat compactions profiles which are closer to those resulting 

from p elastlc than p dlsorder ^ i n f ac t ; numerical simulations have so far described 
with great care the deterministic contact forces (hence the elastic stress) while 
they completely discarded the stochastic processes associated with the explo- 
ration of the many microstates with contacts but zero contact forces. How to 
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Figure 1: Schematic representation of the clastic pressure (dashed line) and dis- 
order pressure (plain line) as a function of the compaction in the range between 
the random loose packing 4> m and the random close packing <t5 m . The disor- 
der pressure vanishes below <p m and diverges at 4>m while the elastic pressure 
vanishes below <f> c and strongly increases above. 



introduce these stochastic processes with the minimum change for the equations 
of motion ? A possibility is to write the equation of motion of particle a in the 
form 

™ a ^3~ = 2 F a(i + F" + m *9 ( 13 ) 

where F af3 is the force exerted at contact with particule /? while the stochastic 
force F a plays the role of the brownian force for thermal systems and witnesses 
to the random exploration of microstates. When tapping granular matter one 
creates a series of short-lasting stochastic forces, while in steady flows this force 
is steadily acting. There is no theory of F a we are aware of but, much like 
the brownian diffusion force is the coarse-grained consequence of the stochastic 
brownian force, the gradient of the disorder pressure acts as a diffusion force 
which is the coarse-grained consequence of the granular stochastic forces. Hence, 
the concept of a rate-independent granular pressure can be used without mod- 
eration in the dense regime, provided this use is restricted to steadily flowing 
media. 
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Figure 2: Dependence of h s top/D on the inclination 9. Comparison between 
experimental results (points) for sand over carpets (from jlfij ) and fitting curves 
deduced from the granular pressure (fl*2"|) (plain) and lfTT|l (dashed) . Experimen- 
tal results with h stop /D < 3.6 were discarded and the value 6 max w 36.1 was 
adopted. 
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